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Abstract 

The simple method for simulation of "thermal" hadron spectra in ultrarelativistic 
heavy ion collisions including longitudinal, transverse and elliptic flow is developed. The 
model is realized as fast Monte-Carlo event generator. 



1 Introduction 



The experimental investigation of ultra-relativistic nuclear collisions offers a unique possibility 
of studying the properties of strongly interacting matter at high energy density. In that regime, 
hadronic matter is expected to become deconfined, and a gas of asymptotically free quarks and 
gluons is formed, the so-called quark-gluon plasma (QGP), in which the colour interactions 
between partons are screened owing to collective effects [1]. One of the important tools to 
study QGP properties is transverse and elliptic flow observables. 

In particular, the experimentally observed growth of the mean transverse momentum with 
increasing mass in most central nuclear collisions at SPS [2, 3] and RHIC [4, 5] energies is natu- 
rally and simply explained with the hydrodynamical model [6] , where the change of momentum 
Apt(i") of a hadron of mass m due to transverse motion of a fluid element at the point r can 
be written as Aprir) = ■msmh.YT{r), where Yt is the collective transverse rapidity. 

Moreover, a strong interest in azimuthal correlation measurements in ultrarelativistic heavy 
ion collisions has recently gained impetus. Recent anisotropic flow data from SPS [7] and 
RHIC [8, 9, 10] can be described well by hydrodynamical models for semi-central collisions and 
transverse momentum, px, up to ~ 2 GeV/c (the coefficient of elliptic flow V2, which is deflned 
as the coefficient of the second harmonic of the particle azimuthal distribution with respect to 
the reaction plane, appears to be monotonously growing with increasing px [11] in this case). 
On the other hand, the majority of microscopical Montc-Carlo models underestimate the flow 
effects (see however [12]). Measurements of V2 present one of the effective tools to test various 
models. In particular, the saturation and gradual decrease of V2 for Pt > 2 GeV/c, predicted 
as a signature of strong partonic energy loss in a QGP, is supported by the recent RHIC 
data [8, 9] extending up to px — 10 GeV/c. The interpolation between the low-px relativistic 
hydrodynamics region and the high-p^ pQCD-computable region was evaluated in [13]. 

The initial gluon densities in Pb— Pb reactions at ^snn = 5.5 TeV at the Large Hadron 
Collider (LHC) are expected to be signiflcantly higher than at RHIC, implying much stronger 
QGP effects. The probing experimental capabilities of LHC detectors together with physics and 
software validation of various Monte-Carlo tools and cross-comparisons among different codes 
are important tasks in the light of coming LHC Heavy Ion Program [14]. A number of Monte- 
Carlo generators is available at the moment to generate heavy ion events at LHC energies: 
HIJING [15], FRITIOF [16], LUCIAE [17], DPMJET-III [18], PSM [19], NEXUS [20], etc. 
However, ffow effects in almost of such models are lacking or implemented insufficiently. Besides, 
running these codes at LHC energies consumes much computing efforts. On the other hand, 
macroscopic hydrodynamical models basically reproduce the bulk of hadron spectra observed 
at SPS and RHIC, and can be in principle used to estimation of particle flow effects at LHC, 
may be extending this approach to even some higher p^ values. Of course, for more detailed 
simulation, one has to take into account the interplay between hydro flow and semi-hard particle 
flow due to parton energy loss, secondary scatterings, etc. 
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2 "Thermal" model and fast Monte-Carlo generation 



Wc suggest simple hydrodynamical Monte-Carlo code [21, 22] giving final hadron spectrum as 
a superposition of a thermal distribution and a collective flow [23, 24, 25, 11], 

E-J^= / f{x,p)p^da, , (1) 



where the invariant distribution function f{x,p) is taken in the Bose- Einstein form for particles 
of integer spin and in the Fermi- Dirac form for particles of half-integer spin (p^ is the 4- 
momentum of hadron, and E — pois its energy). Integration is performed over the hypersurface 
a at the "freeze-out" temperature T — Tf. The formation of the cylindrically symmetric hot 
matter expanding preferably along the cylinder axis is expected in the case of relativistic heavy 
ion collisions; as to the transverse motion, it can be taken into account as a correction [26] . In 
this case, the variables r, r, rj and $ (r = y/x'^ + y^, r = \/t^ — z^, = | In l^^, tan ^ — y/x) 
are commonly used instead of the Cartesian coordinates x, y, z. We consider charged and 
neutral pions, kaons and nucleons only, and kaons and nucleons are supposed to be thermally 
suppressed by their heavier mass. In addition, the linear transverse velocity profile specification 

• 

u — smh Yt = ^— t; 

dr R 

(which follows from a solution of the nonrclativistic continuity equation with uniform density) 
and longitudinal velocity specification in accordance with one-dimensional scaling solution Yi = 
Tj are assumed (Yt and Y^ arc, respectively, the transverse and the longitudinal rapidity of 
collective motion, while R is the effective transverse radius of the system). 

The following procedure were applied to simulate "thermal" hadron spectra in heavy ion 
AA collisions at given impact parameter h. 

1. The 4-momentum p* of a hadron of mass m was generated at random in the rest frame of a 
liquid element in accordance with the isotropic Boltzmann distribution 

f{E*) oc E*y/E*-^ -w?e^^ {-E*/Tf), - 1< cos^ < 1, < 0* < 27r , (2) 

where E* — Vp*^~+rr? is the energy of the hadron, and the polar angle 9* and the azimuthal 
angle 0* specify the direction of its motion in the rest frame of the liquid element. 

2. The spatial position of a liquid element and its local 4- velocity were generated at random 
in accordance with phase space and the character of motion of the fluid: 



/(r) = 2r/R} {0<r <Rf), - r]^^^ < rj < rj^^y,, < $ < 27r, 
sin! 



Ur — sinhYr^^, Ut — Jl + u^ cosh 77, u^ — Jl + sinh 77 , (3) 
Rt 



where -R/ is the effective final transverse radius of the system, which is fixed here by specifying 
the mean charged multiplicity per unit rapidity interval, < dN/dy^ >, in the final state; 
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Vmax — Yf^^ and are maximum longitudinal and transverse collective rapidities. 

3. Further, boost of the hadron 4-momentum in the cm. frame of the event was performed: 



Px 



p* sin 9* cos (f)* + Ur cos $ 



Py — p* sin 6* sin 0* + sin $ 



E* + 
E* + 



Pz — P* COS 9* + Uz E* + 
E = E*Ut + {uY'), 



Ut + 1 



where 



{u^p*^) = UrP* sin 6'* cos ($ - 0*) + u^p* cos 6'* . 



(4) 



(5) 



Anisotropic flow is introduced here under simple assumption that the spatial ellipticity of 
"freeze-out" region, e—{y^ — x^) / (y^ + x^), is directly related to the eUipticity of the system 
formed in the region of the initial overlap of nuclei, Eq — b/2RA {Ra is nucleus radius). This 
"scaling" enables one to avoid introducing additional parameters and, at the same time, leads 
to an azimuthal anisotropy of generated particles due to dependence of effective final radius 
Rf{b) on the angle $ [27]: 



= = 0) min{y 1 - eg sin^ $ + eo cos$, y 1 - eg sin^ $ - eo cos $}. (6) 

Obtained in such a way azimuthal distribution of particles is described well by the elliptic form 
for the domain of reasonable impact parameter values. 

We also set the Poisson multiplicity distribution and assume that the mean multiplicity of 
particles is proportional to the nuclear overlap function [27]. For estimated "freeze-out" pa- 
rameters - temperature Tf — 140 MeV, collective longitudinal rapidity y^"^ — 5 and collective 
transverse rapidity Y^"^ = 1 - we get average hadron transverse momentum < pj^ >— 0.55 
GeV/c and following particle ratios: 

: = 1:1, 



TT^ : : p^ = 24 : 6 : 1, 



71 : 71 



2 : 1, 



p : n = 1 : 1 



The model has been reahzed as fast Monte-Carlo event generator, and corresponding Fortran 
routine is available by the web [28]. The following input parameters should be specified by 
user to set hadron event configuration: beam and target nucleus atomic number; type of event 
centrality generation (options "fixed impact parameter" or "impact parameter is generated with 
standard Glauber geometry between minimum and maximum values" are foreseen); baseline 
mean charged particle multiplicity per unit rapidity at mid-rapidity, < dN^/dy'^ > {y^ = 
0), in central Pb-Pb collisions (total multiplicity for other centralities and atomic numbers 
will be calculated automatically). Since the output particle information is stored in common 
block LUJETS of JETSET routine [29], main users program should be compiled with JETSET 
Fortran routine with extended size (up to 150000) of LUJETS arrays. 
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3 Conclusions 



The simple model to simulate flow effects in heavy ion collisions at LHC energies has been 
developed. This model is realized as fast Monte-Carlo event generator, and corresponding For- 
tran routine is available by the web. 

To conclude, let us to discuss the physics validity of the model application. 

• Internal parameters of the routine for flow were selected as an estimation for LHC heavy 
ion beam energies. The result for other beam energy ranges, obtained without additional 
internal parameters adjusting, is not expected to be reasonable. 

• Hydro-type description of heavy ion collisions is expected to be applicable for central and 
semi-central collisions. The result obtained for very peripheral collisions (6 ~ 2i?^) can 
be not adequate. 

• Hydro flow mechanism in heavy ion collisions is valid for restricted kinematic range: mid- 
rapidity, low and intermediate px- The model is not applicable for very forward rapidity 
{\y\ > 3) and very high (> 2 - 5 GeV/c). 
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